Dark Fluid: a complex scalar field to unify dark energy and dark matter 
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In this article, we examine a model which proposes a common explanation for the presence of 
additional attractive gravitational effects - generally considered to be due to dark matter - in 
galaxies and in clusters, and for the presence of a repulsive effect at cosmological scales - generally 
taken as an indication of the presence of dark energy. We therefore consider the behavior of a 
so-called dark fluid based on a complex scalar field with a conserved t/(l)-charge and associated to 
a specific potential, and show that it can at the same time account for dark matter in galaxies and 
O i 1 in clusters, and agree with the cosmological observations and constraints on dark energy and dark 

matter. 
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I. INTRODUCTION 

After many years of considering that the Universe is filled with a dark matter made of Weakly Interacting Massive 
Particles (WIMPs) and with a dark energy causing an acceleration of the expansion of the Universe, we have still no 
direct evidence of their existences. Moreover, observations of supernovae of type la tend to cause trouble to usual 
dark energy models 1] , and open the way to new kind of models or analyses to explain the observed repulsing effects 
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Alternative models exist to explain the cosmological observations, and in particular some of them try to solve the dark 
energy and dark matter problems by unifying both components into a single "dark fluid" . We can for example note 
that the Generalized Chaplygin Gas model follows this idea and is presently under scrutiny. We have shown in [j| 
that building such a dark fluid model is very difficult, as the model has to be in agreement with many observational 
constraints, and especially has to explain at the same time the cosmological repulsive effects and the local binding 
gravitational effects in the recent Universe. As scalar field-based models for dark energy and dark matter exist in the 
literature Irl M. Isl IqL Eot . it seems interesting to study unifying dark fluid models based on scalar fields. This idea has 
been proposed in The crucial question however concerns the form of the potential of the scalar field. 
In this article, we will consider a complex scalar field and propose an adequate form for its potential. We will show 
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that this scalar field can potentially explain correctly the observations of galaxy rotation curves of spiral galaxies and 
the presence of strong binding gravitational effects in clusters. We will also consider the cosmological behavior of the 
dark fluid and show that it agrees with the cosmological constraints and observations. We will finally conclude by 
suggesting some further directions of investigation beyond this study. 

II. DESCRIPTION OF THE MODEL 

We consider here a scenario in which the dark fluid can be accounted by a complex scalar field associated to a 
conserved charge. We focus on the idea that the scalar field evolves with a quasi-homogeneous density in the Early 
Universe, then produces cluster and galactic halos through Bose condensation and provides repulsing effects outside 
the Bose condensates. The Lagrangian density of the scalar field reads 

C=g^d^d v <j ) -V(4>) , (1) 

and we assume that its potential V is invariant under the global symmetry 

$ -> $' = e l9 $ . (2) 

The principal issue to elaborate a dark fluid model is the choice of a potential. Considering ||, we can conceive that a 
model for a dark fluid would be promissing if it has a negligible density at the Big-Bang Nucleosynthesis (BBN) time, 
if it provides a matter behavior at the time of recombination and of structure formation, and if today its behavior is 
cosmologically repulsive, and attractive in clusters and galaxies. This kind of behavior seems however very difficult 
to achieve. 

We can nevertheless consider the studies of complex scalar fields with quadratic and quartic potentials which were 
performed in 0E1G3. 

In those articles, we showed that these potentials leads to a matter behavior from the time 
of the recombination until today, and are in agreement with the constraints of BBN ^5|. Hence, we can consider 
a potential containing a quadratic term, which will provide a matter behavior. However, this potential has to be 
modified to give a repulsive effect today. Several quintessence potentials, such as a decreasing exponential [16^ or an 
inverse square law |17| . could provide an adequate behavior. In particular, the decreasing exponential seems adequate, 
and the potential reads in this case 

V((f>) = m 2 0t0 + aex p(-/3|0|) . (3) 

a and (3 can then be chosen, so that the quadratic term is dominant until the late time of the structure formation, and 
the exponential term becomes non-negligible only later, once <f> has become small. In that way, it seems possible to 
have a correct cosmological behavior for the fluid. In galaxies, the densities are quite large, and so is <f>, consequently 
the exponential term can again be negligible. In the interstellar medium, where densities are small, </> is also small 
and the exponential term can dominate in the potential and provide a repulsive behavior. 

This potential could then seem interesting to study, but problems could appear when trying to explain the attractive 
gravitational effects at the same time in galaxies and in clusters, due to the different scales involved. Another 
question would concern the way to constraint the parameter /3, which is, apart from the fact that it has to be small, 
quite unclear. Another choice of potential may then be more appropriate. 

Thus, we will rather consider the following potential, constituted of a quadratic term and a Gaussian term: 

V{<t>) = m 2 tf<f> + aexp(-0<f>) . (4) 

This potential can be motivated in the scope of the non-perturbative renormalization group [Tsl ]. The interesting 
feature of this potential is that the second term has the same kind of behavior than the usual decreasing exponential 
potential, but we now have a way to determine [3. In fact, when cf> is small, for example outsi de g alaxies, this potential 
can be approximated by a quartic potential, which is known to lead to a matter behavior [l4|. Therefore, a way to 
constraint (3 would consist in choosing it so that the approximate quartic term leads to a matter behavior at cluster 
scale, whereas the quadratic term leads to a matter behavior at galaxy scale. 

Let us now review the constraints on the three parameters to, a and (3. We will use for this study the re- 
sults of @, m in . 
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Galactic scale: determination of m 

First, we want the mass parameter m to give a matter behavior in galaxies. The Compton wavelength 
associated to the quadratic term is 

H 

'compton = ~ • (5) 



At local scales, we assume that the complex scalar field has a conserved ?7(l)-charge (see Section HTT1 for more 
details) and takes the form 

(7 \ X ) 

</>(&, t) = exp(-iujt) , (6) 

where a is the norm of the scalar field, u> is a constant phase, t is the time and x the position. With a quadratic 
potential, it was shown in 12] that the size of the bosonic halo is then given by 



1 ~ V GVVo Zcom P ton _ V GWq mc ' (7) 
where <jq is the value of a at the galactic center. The typical orbital velocity in such a halo is 

- ~ VG^i . (8) 
c 

For a typical velocity v of the order of 100 km/s and a typical halo size I of 10 kpc, we can estimate the mass 
to be about to ~ 10~ 23 eV, which is in agreement with the results of [l^ . 



• Cosmological scales: determination of a 

We want that today the exponential term of the potential dominates, so that a repulsive behavior could 
dominate at cosmological scales. From the study of the dark fluid parameters performed in 0, it is convenient 
to choose 

a = ^dark energy = °' 71 x Po - 7 x 10" 27 kg m -» , (9) 
where p$ is the cosmological critical density, so that the exponential term plays the role of dark energy. 



• Cluster scale: determination of (3 

Finally, we would like that the exponential term also gives a matter behavior in clusters, i.e. acts as a 
quartic potential when the field has a small norm. In this case, the potential reads, at fourth order in <j>; 

V{<t>) w m 2 t + a-a/30 t ^+^a/3 2 (0 t 0) 2 (10) 
m a + (m 2 - m' 2 )^ ' <f> + A(0 t 4>f , (11) 

with ml 2 — a/3 and A = a(3 2 /2. We can assume here that m > to'. We have shown in 0] that the scalar field 
condenses on a distance L, which we would like to be of the typical size of a cluster (L»l Mpc). In this case 

8TTG?7l 4 rrL 2 

A = f£— , (12) 

where w 2 g- = m — to' 2 . The typical cluster scale of 1 Mpc then imposes 

A-10" 89 , (13) 
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and 

0= „ eV- 2 . (14) 

V a 

We can then calculate m' to confirm that m > m', and show that the effective mass is not very different from 
the quadratic mass (m e g f» 0.999 m > 0). 



Thus, with these choices of parameters, the potential has a single non-degenerate minimum for <fi = 0, which eludes 
the question of global cosmic strings. 

Considering that the different parameters of the model have been fixed, we can now analyze the behavior of the field 
at galactic, cluster and finally cosmological scales for these parameters. 



III. GALACTIC SCALE 



We will first show that this scalar field can play the role of dark matter at galactic scale. Let us consider that a galaxy 
evolves in a space represented by a static and isotropic metric 

dr 2 = e 2u dt 2 - e 2v (dr 2 + r 2 d9 2 + r 2 sin 2 9dcj> 2 ) , (15) 

where u and v only depend on the position. We will assume a spherical symmetry for the system, as it was already 
shown that this approximation leads to correct and significant results [TJ, and study the rotation curves generated 
by the model for spiral galaxies. 

The scalar field is as usual to be considered to have an internal rotation, so that we can write, in a spherical 
symmetry, (f>(x,t) cx a(r 1 t)e~ ie ^' r ' t \ where 9 is the internal rotation angle. However, a deeper study of the cosmolog- 
ical behavior of 9 (see in particular equation (|42|) . in which u — 9) could reveal that when the scalar field behaves 
like cosmological matter (i.e. when a 2 oc a -3 , with a the cosmological scale factor), the first derivative of 9 with 
respect to time (u>) is nearly a constant. We thus assume that <f>(x,t) ~ <r(r, £)e -lw(r - ) *. As we study a static model, 
we should also consider that the modulus of the scalar field is time-independent. Besides, we will first assume that 
u>(r) is space-independent (but galaxy-dependant). We will discuss this approximation at the end of this section. 
Thus, following these assumptions, the scalar field has a uniformly rotating phase u>, which corresponds to the 
simplest possible realization for a scalar field. One could refer to ^} for an alternative analysis of the behavior of a 
similar scalar field in galaxies. We can write 

<t>{x,t) = ^-exp(-iwt) . (16) 
v2 

In the case where the internal rotation period 2tt/lo is much smaller than the rotation period of the galaxy, the 
Klein-Gordon and Einstein equations read 

e~ 2v | a" + (u' + v' + ^ a' | c 2 + oj 2 e - 2u a - V'(a) = , (17) 

2v" + v' 2 + = —e 2v lp b + e 2u — — + e 2v — + V \ (18) 

and 

II , II , 12 , 1 / / . I\ 8n G 2v \ -2u Uj2(t2 -2v a ' 2 Tr\ rin s 

u +v + u +~(u + v) = — g-e < e — e — V > , (19) 



where 



1 22 , „ ( 1 a 2 



V(a) = 2™ CT +«exp ) (20) 

V'(a) = m 2 a - a(3aexp(-^f3a 2 ) . (21) 
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We can recover the Newtonian limit by posing u — $ and v = — ty, where <E> is the usual gravitational potential. A 
careful examination of a' [l4| can reveal that a' can be safely neglected in comparison to loo . 
In this limit, the scalar field contributes to the mass of the galaxy by adding an effective density 

p eS = 4^0 - 2V(<p) = 2wV - m 2 a 2 - 2acxp (^-^P^ 2 \ (22) 

to the baryon density p b . The effective pressure of the scalar field reads 

P e& = g ^^-V^) = \^-\m^-ae^(-\pA . (23) 

In the Newtonian limit, <I> = "J, so that equation (|I9(I leads to 

uj 2 <t 2 w m 2 a 2 + 2aexp ■ (24) 

Thus, the effective pressure vanishes, and the scalar field behaves like pressureless matter. 

In spiral galaxies, the typical density at large radii is of the order of 10~ 21 kg m~ 3 , which corresponds to 
a ~ 10 -2 kg 1 / 2 m~ 3 / 2 s. Therefore, in this case, a exp (— \fio 2 ^) ~ 10~ 3000 kg m -3 is completely negligible, and 
only the quadratic term remains. Thus, we are in the same case as the analysis of [12J , which considered a simple 
quadratic potential. We have then uJ 2 m m 2 and consequently 

p cff «mV . (25) 

For m = 1CP 23 GeV, we can calculate the internal rotation period, and show that it is of the order of 30 years, which 
confirms the hypothesis that the internal rotation is much faster than the galactic rotation. Concerning the baryons, 
should the disk be alone, the rotation curve would be 

a, \ 2 ( v l-97(r/r opt ) 1 - 22 

{{r/r opt y + 0.78^} 1 - 4ci 

leading to a baryonic density 

, , 1 v 2 d (r op t) ( 4.38(r/r opt )^- 78 5.64(r/r opt ) 122 \ 

Pb[ 1 AttG r 2 opt \{(r/r opt ) 2 + 0.78 2 }^ {(r/r opt ) 2 + 0.78 2 } 2 « J ' ^ ' > 

r pt being the optical radius, defined as the radius of the sphere encompassing 83% of the luminous matter. 



When solving the equations of motion, we can obtain the result that, for each set of parameters, the system 
can have a discrete number of possible configurations, corresponding to fundamental and excited state of the Bose 
condensate [l2T |. We can derive the rotation curves from the Newtonian gravitational potential. 

Restricting us to the fundamental state, it is possible to generate approximately flat rotation curves u p to large radii, 
as shown on figure ^ As the potential in galaxies restricts to the quadratic term, we can refer to \\2\ for a more 
complete study of such a potential at galactic scales. 

The size and the density of the halos are fixed by the value of ui. Then, to each galaxy corresponds a specific u>. As 
the Newtonian approximation leads to lu 2 m to 2 , to can only slightly vary, which justify the previous approximation 
that lu is space-independent inside a galaxy. However, if we want to study the variation of this u> between galaxies, 
we have to consider an extra equation coming from the imaginary part of the Klein-Gordon equation: 

J'ar + 2u'a'r + {2 + r(u' + v')}oj'<j = . (28) 

We can nevertheless note that, because of the very small variations of cj, this equation has no effect inside a specific 
galaxy, but can however explain the presence of different effective values of uj inside a galaxy cluster. 

To conclude this section, this scalar field may account for dark matter in galaxies, what is confirmed in |lOj |. 
Let us now consider the behavior of the scalar field at cluster scale. 
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FIG. 1: Rotation curve of a generic spiral galaxy, in presence of a complex scalar field with a uniformly rotation phase 
accounting for dark matter. The solid line corresponds to an example of rotation curve provided by the complex scalar field 
plus the baryons, whereas the dashed line shows the contribution of the baryons only. 

IV. CLUSTER SCALE 

In clusters, one can expect the same kind of behavior as in galaxies, an important difference being that the average 
densities are smaller. We can consider that the typical density in clusters is around 1CP 25 kg m~ 3 . From equation 
B22JI . we can show that in the cluster case, cvexp (— ^/3er 2 ) ~ a. This term cannot be neglected anymore, but can be 
approximated by a — \afio 2 + ^a(3 2 a A . The potential then reads: 



V{a) wa + ^(m 2 - m! 2 )a 2 + ^Xa* = a+ ^m| ff a 2 + ^Xa 4 



where m' 2 = a(3, m 2 



off 



ml 2 and A = af3 2 /2. In this case, the Klein-Gordon equation becomes 



Xa 2 = (1 - 2$V 2 - m 



off 



(29) 



(30) 



This result is similar to the one derived in [l4| . We can remark that the Klein-Gordon equation does not depend on 
a. The effective density of the scalar field rewrites 



Pcft 



o 2 2 2 2 

= Iuj a — m eff <7 



T Xa 4 



(31) 



The scalar field is then in a Bose condensate state, on the boundaries of which the gravitational potential is 

1 



$n = - 1 - -ir 



(32) 



As the potential <&o is small, we have lu 2 rs m^g, and the effective density becomes 

2m 4 



Pcff 



-2a 



X 



-($ -$)W($o-$) 



The Klein-Gordon and Einstein equations can then be combined into the following Poisson equation |l4 

A$ = 4irG( Pb + p cS ) . 



(33) 



(34) 



Hence, the scalar field has a matter behavior up to the end of the Bose condensate. The value of A has to be chosen so 
that the size of the condensate is approximately the size of a typical cluster, in the way presented in the introduction. 
In the case where the luminous matter density can be neglected in comparison to the density of the scalar field - 
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which is a decent hypothesis at cluster scale and if the halo can be considered spherical, the Poison equation reduces 
to the Lanc-Emdcn equation for a polytrop n=l. We can therefore determine analytically the effective density of the 
scalar field: 

p = -2a + p° — , (35) 

where p° is the scalar field density at the center of the cluster, and z = r/L, L being the typical size of a cluster. The 
effective size of the bosonic halo is then: 

r \ 1 1/2 1 

*=*kc} 55- (36> 

and the total mass of the cluster is determined by the value of a at the center. 

Knowing that p^ « 1CP 25 kg m~ 3 and a s=a 7 x 10~ 27 kg m~ 3 , the exponential term will have only a slight effect inside 
the Bose condensate, but far away from the center of clusters, beyond the end of the Bose condensate, only this term 
remains in the effective density: 

P eS « - 2 " • (37) 

Consequently, outside clusters, the behavior of the dark fluid is repulsive, and so completely similar to that of a 
cosmological constant. 

In summary, we can interpret our results in this way: where baryons are gathered in sufficient densities, the 
scalar field is excited and condenses into attractive halos, and therefore has a matter behavior, whereas for a low 
baryonic density, the field is not excited anymore and has a repulsive behavior. We can thus guess that the global 
behavior of the field is repulsive, in particular at cosmological scale. We will verify this idea in the next section. 



V. COSMOLOGICAL BEHAVIOR 



Let us now consider the cosmological behavior of the scalar field, in a flat, homogeneous and isotropic Universe, filled 
only with radiation, baryonic matter and the dark fluid scalar field. The pressure and density of the scalar fluid can 
be written 



and 



where the scalar field has an internal rotation 



^ct> + v(4>) , 



The Klein-Gordon and Friedman equations read 



d 2 a 3 da da 2 
dt A a dt dt 



(3a 2 - u?a = 



duj 3 da da 

—a H — uoa + 2uu— = 

dt a dt dt 



(38) 
(39) 

(40) 

(41) 
(42) 



3H A 



da 
a dt 



8ttG 



Prad + Pb + 



2 2, 2 2 I 

+ lo a + m a > + a exp 



-Ba- 



rn 



where p rac i is the usual density of relativistic photons and neutrinos, p^ is the usual baryonic density, and u> = d6 /dt 
is time-dependant. The second equation implies the conservation of the charge per comoving volume Q = uia 2 a 3 . 
Therefore, we can rewrite the first equation as 



d 2 a 3 da da 



-aW + alTtTt +m(J - a(3(JC ^\-^ I- 



1 



Q 2 



= o . 



(44) 
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Using the same parameters as defined before, we can solve the equations and determine the evolution of the cosmo- 
logical parameters Q = p/p$ in function of the scale factor a. The initial value of the field density can be chosen 
arbitrarily and leads to an identical global behavior. The charge Q has been fixed so that the equality scalar field 
density/radiation density occurs approximately at the same time as the equality dark matter/radiation in the stan- 
dard model of cosmology. The cosmological evolution of the density of the scalar field in this context is plotted on 
figure El 

One can remark that the field follows a four-step evolution, whose first three steps had already been observed in [l3|. 




FIG. 2: Cosmological evolution of the density of the scalar field dark fluid in comparison to the densities of baryonic matter 
and radiation. 



and which is directly related to the form of the Klein-Gordon and Friedman equations: 

• The first behavior, at very high redshift, occurs when the kinetic term of the scalar field dominates, leading to 
a a~ 6 decay of the scalar field density. At this time, the scalar field density is very high, and the potential has 
a negligible influence on the evolution of the field. This step is in fact determined by the high expansion rate of 
the Universe, so that only the first two terms of equation l|41|) dominates. As at the BBN time the field density 
is negligible in comparison to the radiation density, this behavior is in agreement with the BBN constraints. 

• This behavior stops as the expansion rate decreases, i.e. as the influence of the potential increases. A transition 
will occur when the potential will exactly compensate the kinetic terms in equation (|41|l . The scalar field density 
will then reach a plateau, leading to a cosmological constant behavior. 

• This cosmological constant behavior will continue until H ~ m. At this time, the mass term becomes important 
in the potential, and the field oscillates quickly in time. Its dynamics is then comparable to that of an effective 
real scalar field oscillating rapidly, and its time-averaged pressure vanishes. Thus, the field then behaves like 
pressureless matter. This transition occurs before matter-radiation equality. 

• Recently, at a late time in structure formation, as the densities continue to decrease and the exponential part of 
the potential becomes non-negligible, the oscillations have damped sufficiently so that the field can settle into 
its non-degenerate vacuum. This leads to a dark energy behavior at low redshifts. However, the transition to 
this dark energy behavior is not completed yet, and in the future, one can expect the scalar field to achieve its 
transition and behave like a cosmological constant. 

This scenario is in complete agreement with the constraints on a unifying dark fluid model which were derived in |j| . 



VI. CONCLUSION 



In this work, we have shown that it is possible to elaborate a dark fluid model replacing the standard dark energy / 
dark matter models, using a complex scalar field. We have seen that a potential containing a quadratic term and a 
Gaussian term is suitable to explain the gravitation observations at galaxy, cluster and cosmological scales. Based on 
|12|. we also know that this model can account for a large variety of rotation curve shapes. Thus, we can conclude 
that this model could be an interesting alternative to the dark matter / dark energy models. 
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Of course, many studies are still needed to test this model. In particular, it would be interesting to consider the 
formation of large scale structures and to study the influence of the dark fluid on the anisotropies of the Cosmic 
Microwave Background. Another question concerns the form of the potential. Indeed, this potential has interesting 
properties, but is not directly derived from high energy theories. We can consider it as a toy-potential rather than a 
definitive answer. It would then be interesting to perform the study of different kinds of potentials, and perhaps to 
find an adequate potential derived from a high energy theory. 

In any case, this simple model looks promissing, and reinforces the suggestion that dark fluid models should be studied 
deeply. 
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